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We derive a coordinate-independent formulation of the post-l-Newtonian approximation to gen- 
eral relativity. This formulation is a generalization of the Newton-Cartan geometric formulation of 
Newtonian gravity. It involves several fields and a connection, but no spacetime metric at the fun- 
damental level. We show that the usual coordinate-dependent equations of post-Newtonian gravity 
are recovered when one specializes to asymptotically flat spacetimes and to appropriate classes of 
coordinates. 
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I. INTRODUCTION AND SUMMARY 
A. Background and Motivation 

The weak field, slow motion approximation to general 
relativity, also called the post-Newtonian approximation, 
consists of expanding in the small parameters v 2 /c 2 and 
$/c 2 , where v is a typical velocity of the system un- 
der consideration, $ is the Newtonian potential, and c is 
the speed of light. At the leading order Newton's the- 
ory is recovered, and higher order corrections are called 
post-l-Newtonian corrections, post-2-Newtonian correc- 
tions and so on. This approximation scheme is very useful 
in astrophysical applications and is very well developed. 
Reviews can be found in Ref. [l[ and in the book by Will 
@, and a historical review can be found in Ref. 

There are two types of equations that arise in post- 
Newtonian theory. The first are continuum equations of 
motion, for example for gravity coupled to a perfect fluid, 
for which one obtains generalizations of the equations of 
Newtonian hydrodynamics. Such continuum equations 
have been used extensively in numerical simulations (al- 
though fully relativistic simulations are now the state of 
the art The second type of equations are "point 

particle equations" , which describe the motions of bod- 
ies whose sizes are small compared to their separations. 
Such point particle equations (and extensions to include 
spins) can be derived from the underlying continuum the- 
ory by a variety of methods. Currently the equations of 
motion for two point particles are known up to post-3. 5- 
Newtonian order, see, for example Ref. Q and references 
therein. In this paper we shall be concerned only with 
continuum equations. 

Over the years, the foundations of the Newtonian and 
post-Newtonian approximations have been studied in de- 
tail and with considerable mathematical rigor by a num- 
ber of researchers. Futamase and Schutz [6] have shown 
that the various orders of post-Newtonian approxima- 
tion are asymptotic approximations to fully relativistic 
solutions. Frittelli and Reula [3] showed that, given a 
solution of the equations of Newtonian gravity, there ex- 



ists a one-parameter family of exact relativistic solutions 
which for a finite amount of time are close to the Newto- 
nian solution. Rendall Q gave a mathematically rigorous 
derivation of both the Newtonian and post-Newtonian 
approximations from a precise set of axioms. 

However, there remains one aspect of post-Newtonian 
theory which has not been fully explored: there is as 
yet no covariant version of the theory. Usually, in or- 
der to find the equations of post-Newtonian theory, one 
first introduces specializations of the coordinate system 
(or gauge), and then Einstein's equations are expanded 
order by order. The gauge specializations are chosen to 
simplify the resulting equations. For example, at post- 
l-Newtonian order, the harmonic gauge condition and 
the so-called standard post-Newtonian gauge condition 
Q are often used. In each of these gauges the post- 
Newtonian equations take a different form. The situa- 
tion is similar to knowing the laws of electromagnetism 
only in a handful of gauges (e.g. the Lorentz and the 
Coulomb gauge), without knowing the underlying gauge 
independent equations. 

This lack of covariance of post-Newtonian theory as 
currently formulated has several disadvantages: 

• If one is attempting to compare two different calcu- 
lations, it is often helpful to identify and compute 
gauge-invariant quantities. Generally, such quanti- 
ties are easier to identify starting from a covariant 
formulation of the theory. 

• In attempting to developing intuition from the 
post-Newtonian equations, it can be difficult to sort 
out which aspects of the equations contain the ac- 
tual gauge-independent physics and which aspects 
are gauge dependent. For example, there is a well- 
known analogy between the equations of post-l- 
Newtonian theory in certain gauges and those of 
electromagnetism. In this case the analogy with 
electrostatics and magnetostatics is physical, but 
the additional aspects of the analogy concerning 
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magnetodynamics are gauge. 1 

• At post-l-Newtonian order, there is a well- 
developed framework for celestial mechanics that 
describes N interacting, deformable bodies 
This framework is quite complicated, involving sep- 
arate coordinate systems for each body as well as 
a global coordinate system, together with a set of 
fields associated with each of the coordinate sys- 
tems. A covariant formulation might simplify some 
aspects of this framework. 

B. Covariant post-Newtonian theory 

The purpose of this paper is to derive a coordinate in- 
dependent formulation of post-l-Newtonian theory. We 
will derive a fully covariant set of equations, involving a 
number of tensor fields and a connection, which reduce to 
the standard post-Newtonian equations in specific coor- 
dinate systems for asymptotically flat spacetimes. In the 
case of Newton's theory, such a geometric formulation 
has already been found by Cartan and others [l3l - [r7j ]. 
and is called Newton-Cartan theory. Building on earlier 
work of Dautcourt [H, HU , our derivation will reproduce 
Newton-Cartan theory at leading order, and at the next 
order will give a covariant version of post-l-Newtonian 
theory which we call "post- Newton-Cartan theory" . 

We will actually derive two different versions of post- 
Newton-Cartan theory. The first version, which we call 
perturbative post-Newton-Cartan theory, allows one to 
compute the leading order corrections to a solution of 
Newton-Cartan theory. It is a unique theory, in which 
post-Newtonian corrections to Newtonian quantities are 
treated as independent variables to be solved for. The 
second version, which we call combined post-Newton- 
Cartan theory, combines some of the Newtonian and 
post-Newtonian variables together. It is more econom- 
ical and convenient to use than the first version, because 
it has fewer variables and fewer equations. Solutions of 
this theory will be accurate to post-l-Newtonian order, 
but will in addition contain post-2-Newtonian, post-3- 
Newtonian etc. pieces. The combined theory is not 
unique; different choices could be made to define different 
theories whose solutions differ at post-2-Newtonian and 
higher orders. We note that standard, coordinate-specific 
post-Newtonian theory also comes in "perturbative" and 
"combined" versions |6L 



C. Results 

We now turn to a description of our results. We start 
by reviewing Newton-Cartan theory, then we describe the 



If the analogy with electromagnetism were complete there would 
be radiation in the post-l-Newtonian theory. 



TABLE I: Fields and equations of Newton-Cartan theory. 
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two versions of post-Newton-Cartan theory. 

1. Newton-Cartan theory 

The variables of Newton-Cartan theory are a one-form 
t a , a symmetric contravariant tensor field h ab with sig- 
nature (0,+,+,+), and a torsion-free connection D a . 
Matter is described by a symmetric contravariant stress- 
energy tensor T ab ■ The fields t a and h ab are nondynami- 
cal, background fields, while D a and T ab are dynamical. 

The equations of the theory are the orthogonality con- 
dition 

h ab t b = 0, (1.1) 

the compatibility of the fields h ab and t a with the con- 
nection, 

D a h bc = 0, (1.2a) 
D a t b - 0, (1.2b) 

and the Trautman condition 

h fla R f(bc) d] =°> (°) 

where R abc d is the Riemann tensor associated with the 
connection D a . In addition we have the field equation 

R ab = 4iTt a t b t c t d T cd , (1.4) 

where R ab = R acb c , and the stress-energy conservation 
equation 

D a T ab - 0. (1.5) 

These fields and equations are summarized in Table [I] In 
Sec. lIII Al below we review the derivation of this Newton- 
Cartan theory from general relativity, and in Sec lIVBl wc 
review how Newton-Cartan theory reduces to standard 
Newtonian gravity in appropriate circumstances and in 
appropriate coordinate systems. 
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TABLE II: Additional fields and equations of perturbative post-Newton-Cartan theory. 
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TABLE III: Fields and equations of combined post-Newton-Cartan theory. 



Field 


Description 


h ab 


spatial metric, (0, +, +, +) 
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time one-form 


k ab 


contravariant metric perturbation 
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covariant metric perturbation 
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2. Perturbative post-Newton-Cartan theory 

Turn, now to the perturbative post-Newton-Cartan 
theory. This theory contains the Newtonian fields h ab , 
t a , D a and T ab , and in addition four new fields: a sym- 
metric, contravariant tensor k ab , a symmetric, covariant 
tensor p ab , a perturbation A c ab to the connection, and a 
perturbation S ab to the matter stress-energy tensor. 

The equations of the theory are the six Newtonian 
equations (jl.ip - (11.5[) . together with six post-Newtonian 
equations: (i) the orthogonality condition 

h ab Pbc -k ab t b t c = 8 a c ; (1.6) 

(ii) the compatibility of k ab with the connection, 

D a k bc + A b ad h dc + A c ad h bd = 0; (1.7) 

(iii) the compatibility of p ab with the connection, 

D aPbc + A d ab t d t c + A d ac t b t d = 0; (1.8) 



(iv) the post-Newtonian Trautman condition 

fce[ X( 6c) dl - ^AMc + h e[a D [b A d \ )e = 0; (1.9) 

(v) the field equation 

- D a A b bc + D b A b ac = 4n[t a t b t c t d S cd - 4t c t (aPb)d T cd 
+p cd T cd t a t b + t c t d T cd Pab ] ; 

(1.10) 

and (vi) the stress-energy conservation equation 

D a S ab + A a ac T cb + A b ac T ac = 0. (1.11) 

These fields and equations are summarized in Table UT1 In 
Sec. lIII Al below we review the derivation of this perturba- 
tive post-Newton-Cartan theory from general relativity, 
and in Sec. IIV CI we review how it reduces to standard 
post-l-Newtonian gravity in appropriate circumstances 
and in appropriate coordinate systems. 
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3. Combined post-Newton-Cartan theory 

The combined post-Newton-Cartan theory has aspects 
of both the Newton-Cartan and the perturbative post- 
Newton-Cartan theories. The independent variables are 
a one-form t a , a symmetric contravariant tensor field h ab 
with signature (0,+,+,+), a symmetric, contravariant 
tensor k , and a symmetric, covariant tensor p a b- There 
is also a torsion- free connection D a , which is defined ac- 
curate to post-l-Newtonian order. Matter is described 
by a symmetric contravariant stress-energy tensor T ab , 
which is also defined accurate to post-l-Newtonian or- 
der. The fields t a and h ab are nondynamical, background 
fields, while the remaining fields are dynamical. 

The equations of the theory are: (i) the orthogonality 
conditions 

h ab i b = 0, (1.12a) 
h ab p bc - k ab t b t c = 5 a c - (1.12b) 

(ii) the connection compatibility conditions 

D a {h bc + k bc ) = 0, (1.13a) 
Da{-i b t c +Pbc) = 0; (1.13b) 

(iii) the Trautman condition 

(h ela + k ela )R e{bc) d] (D e ) = 0; (1.14) 

(iv) the field equation 

Rab(D c ) = 4n[tJbi c idT cd - 4t c t (aPb)d f cd 

+Pcdt cd i a i h + icidt cd Pab] ; (1.15) 

and (v) the stress-energy conservation equation 

D a t ab = 0. (1.16) 

These fields and equations are summarized in Table 
IIIII The derivation of this combined post-Newton- 
Cartan theory from the Newton-Cartan and post- 
Newton-Cartan theories is given in Sec. IIIIBI below. 

D. Organization of this paper 

In Sec. [n] we list the assumptions underlying our 
derivation, and discuss the motivation for these assump- 
tions. Section IIIII gives the derivation of the Newton- 
Cartan theory and both versions of the post-Newton- 
Cartan theory, from these assumptions. In Sec. IIVI wc 
show that the equations of perturbative post-Newton- 
Cartan theory reduce to the standard coordinate-specific 
equations of post-Newtonian theory, under certain con- 
ditions and in certain coordinate systems. We conclude 
by summarizing our results and their implications in Sec. 

El 



Some technical and side issues are discussed in the 
Appendices. Appendix [X] derives an alternative form of 
one of our assumptions. Appendix [B] derives the gauge 
transformation properties of the fields of the perturbative 
post-Newton-Cartan theory. Finally Appendix lO special- 
izes the formalism, as an example, to perfect fluids. 



E. Notation 

Throughout we will use the metric and sign conven- 
tions of Misner, Thorne and Wheeler [l7|. Furthermore 
we will use Penrose's abstract index notation 20], with 
indices (a, &,...) from the beginning of the Latin alphabet 
denoting general tensors. When we specialize to particu- 
lar coordinate systems, we will use indices (a, j3, . . .) from 
the Greek alphabet to denote general components of ten- 
sors. We will also use indices (i,j,k,...) from the middle 
of the Latin alphabet to denote spatial components of 
tensors, and an index to denote time components, in 
particular coordinate systems. For example, u a will de- 
note a four-velocity vector, u M will be the components 
of the four- velocity in some coordinate system, and u l 
and u° will be the spatial and time components in this 
coordinate system. We will work with units in which 
G = c = 1, except for some special cases in Sec. Ill Bl 
where factors of c and G are explicitly included. Finally, 
we will use the conventional definitions of the two order 
symbols 0() and oQ. 



II. FOUNDATIONS AND ASSUMPTIONS 

In this section we define and discuss the assumptions 
that we make in order to derive the Newton-Cartan 
and post-Newton-Cartan theories from general relativ- 
ity. At Newtonian order we follow closely the treatments 
of Dautcourt [3 El and of Rendall @. 

A. Assumptions 

The starting point is to assume a one-parameter family 
g a b{e),T ab (e) of exact solutions of Einstein's equations 

G ab [g cd (e)] = SnT ab (s), (2.1) 

where e is a real parameter with < e < £q for some 
£q. Wc assume that the solutions are smooth functions 
of spacetime and of e for e > 0. We do not assume the 
existence of a solution at e = 0. 
The key assumptions we make are: 

1. The contravariant metric g ab (e) can be expanded 
near e = as 

g ab (e) = h ab + ek ab + e 2 j ab + e 3 l ab + o(e 3 ), (2.2) 
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where h ab , k ab , j ab and l ab are e- independent sym- 
metric tensor fields on spacetime. Furthermore h ab 
has signature 

fc«*~(0, +,+,+), (2.3) 

and k ab t a tb is everywhere nonzero, where t a is the 
direction defined by h ab tt, = 0. 

2. The matter stress-energy tensor can be expanded 
near e = as 

T ab (e) = e 2 T ab + e 3 S ab + 0(e 4 ), (2.4) 

where again T ab and S ab are e-independent tensor 
fields on spacetime. 

3. The connection V a associated with the metric 
g a b(s) has a continuous limit D a as e — > 0, 



V a = D a + 0(e). 



B. Motivation and Discussion 



(2.5) 



The Newtonian limit of general relativity is the limit c — > 
oo at fixed G. 

The first assumption, the expansion (|2.2j) of the con- 
travariant components of the metric, can now be moti- 
vated as follows. The spacetime metric should be close 
to the flat, Minkowski metric in the Newtonian limit. In 
system- adapted units where T = M. = C = 1, this is 



ds 2 = 
gives 



-c 2 dt 2 



■ dx + dy + dz . Now identifying e = c 



g» v = diag(-e, 1,1,1) = 0(1) +0(e), 



which satisfies assumption 1 and has signature 
(0,+,+,+) at order 0(e°). Corrections from the New- 
tonian potential do not change this conclusion. 

The second assumption, the expansion (|2.4|) of the 
stress-energy tensor, can be motivated similarly using 
dimensional analysis. In a general system of units Ein- 
stein's equation is 



8ttG 



(2.7) 



We now discuss the motivation for and properties of 
these assumptions. 

First, we note that, since the limiting contravariant 
metric h ab is degenerate, the limit as e — > of the covari- 
ant metric does not exist. Therefore there is no limiting, 
background solution of Einstein's equations at e = in 
this framework, unlike the situation for standard pertur- 
bation theory [|, QJ, QJ]. 

Next, the assumptions are explicitly local and covari- 
ant. In Sec. Mil below we will show that the Newton- 
Cartan and post-Newton-Cartan theories can be derived 
from them in a local and covariant manner. Suppose 
now that one demands that the assumptions apply only 
in a given, finite region of spacetime. Then, the Newton- 
Cart an and post-Newton-Cartan equations of Tables U 
and |TT] will be satisfied in that region. However, as is 
well known, it does not follow that the usual equations of 
Newtonian gravity will be satisfied, since Newton-Cartan 
theory contains more local degrees of freedom than New- 
tonian gravity. The physical reason for this will be dis- 
cussed in Sec. IIVB1 below. In order to obtain Newtonian 
gravity, it is necessary to assume an asymptotically flat 
spacetime and to impose the assumptions (|2.ip - (|2.5|) 
throughout all of spacetime [1, [TU, [l9j]. We will find a 
similar situation for the post-Newton-Cartan theory in 
Sec. IIVI below: it contains more local degrees of freedom 
than standard post-l-Newtonian general relativity, and 
reduces to it only when the assumptions (|2.ip - (|2.5[) 
hold globally in an asymptotically flat spacetime. 

Consider now an isolated physical system that is char- 
acterized by some mass scale M, lengthscale C, and 
timescale T . Then, from Newton's constant of gravi- 
tation G and the speed of light c one can form two di- 
mensionless parameters: 



where from dimensional analysis T u ~ M.C, 3 , T u ~ 
MC- 2 T~ X , and T ij - MC- 1 ^ 2 . In particular, all the 
components of T Q/3 are independent of c to leading order. 
If we now specialize to system-adapted units, all of the 
components of T a/3 are of order unity, and the right-hand 
side of Einstein's equation is of order G/c 4 oc e 2 , since 
we are considering a limit in which G is held fixed 2 . 

Before discussing the third assumption, it is useful to 
consider the gauge freedom present in the formalism. It 
may appear that the formalism so far is explicitly gauge 
invariant, since it is covariant under general diffeomor- 
phisms. However, given a one-parameter family of solu- 
tions g a b(s),T ab (e) of Einstein's equations on a manifold 
M, the gauge freedom consists of a one-parameter fam- 
ily of diffeomorphisms ip e : M — > M, which act on the 
solutions via 



9ab{s)^Ve*9ab{s), T a \e) -> Ve , T a0 (e). (2.8) 

Here ip e * is the pullback mapping on tensor fields that 
is defined in, for example, Appendix C of Ref. [20|. An 
important point is that, since we do not require the exis- 
tence of a solution at e = 0, there is no reason to require 
the diffeomorphisms ip £ to have a well defined limit as 
e — > 0. Thus, there are two subclasses of gauge transfor- 
mations: 

• Transformations which we will call regular, consist- 
ing of smooth one-parameter families of diffeomor- 
phisms which have a smooth limit as e — > 0. Such 



cT 
C ' 



G 



GMT 2 
£ 3 ' 



(2.6) 



2 If we replace the assumption (12.41) with a power series expansion 
of T ab that starts with a term proportional to e" with v ^ 2, we 
would obtain a non-Newtonian limit of general relativity where 
c — > oo with Gc 2 " -4 held fixed. 



6 



families can be parametrized in terms of a fixed, e- 
independent diffeomorphism tp and a set of vector 
fields , £" 2 ^ , . . . via the expansion [2l[ 

<p e = <pooV i(i) {e)oV^ 2) (e 2 )o..., (2.9) 

where for any vector field £, T>^(s) is the diffeomor- 
phism given by moving any point e units along an 
integral curve of £. 

• Transformations which we will call irregular, con- 
sisting of smooth one-parameter families ip £ of dif- 
feomorphisms which have do not have a smooth 
limit as e — > 0. 

Our assumptions are explicitly covariant under regular 
gauge transformations, as will be discussed in more de- 
tail in Appendix |B] below. However, they are not covari- 
ant under irregular gauge transformations. For example, 
consider the prototypical Newtonian-order metric that 
satisfies our assumptions (|2.ip - (|2.5p . namely 

ds 2 = -i [l + 2e<S>{t,x t ) + 0{s 2 )} dt 2 

+ [5ij + 0(e)]dxW, (2.10) 

where $ is the Newtonian potential. Under the irregular 
gauge transformation t — > \fet, x % — > x l , this metric is 
transformed into 

ds 2 = - [l + 2e${v^t,x l ) + 0{e 2 )] dt 2 

+ [S ij + 0(E)]dx i dx j , (2.11) 

which does not satisfy our assumptions (|2.1I) - (|2.5I) . 
Thus, our assumptions do entail a certain limited amount 
of gauge specialization 3 , even though they are covariant 
under transformations of the form (|2.9[) . 

We now turn to a discussion of the third assumption, 
the expansion (|2.5p of the connection V a . First, we re- 
mark that assumptions 1 and 2 are insufficient to char- 
acterize the Newtonian limit. For example, consider a 
one-parameter family of static vacuum spacetimes, of the 
form 

ds 2 = -e 2a{xh ^dt 2 + h tJ (x k ,e)dx l dx : >, (2.12) 



Because of the limited gauge dependence of our assumptions, it 
is not a priori obvious that the fields D a , t a , h , k ab , etc. that 
characterize the Newton-Cartan and post-Newton-Cartan theo- 
ries arc physically unique. More precisely, suppose that we start 
with a one-parameter family of solutions g a i,(£),T a (e) which 
satisfies our assumptions, and is thus characterized by a set of 
limiting fields D a , t a , h ab , k ab etc. Now make a general (possibly 
irregular) gauge transformation, to obtain a new one-parameter 
family <? a (,(£), f ab (e) of solutions. If this new family also satis- 
fies our assumptions, then it will be characterized by a new set 
of fields D a , t a , h ab , k . We conjecture that in this case there 
must exist a regular gauge transformation of the form J2 .9D relat- 
ing the two sets of fields (in the manner described in Appendix 
HP , so that the fields are unique in a physical sense. 



which is smooth in e near e = 0. This one-parameter fam- 
ily of metrics, when written in terms of the coordinates 
t = Jet and a;- 5 , satisfies our assumptions 1 and 2, but is 
not of the type associated with the Newtonian limit. In 
particular components of the Riemann curvature tensor 
will diverge in this example as e — > 0. Therefore some 
additional assumption like assumption 3 is necessary. 4 
Our assumption 3 is actually slightly stronger than is 
necessary: we show in Appendix [A] that, whenever as- 
sumptions 1 and 2 hold in a local region, and the Rie- 
mann tensor R abc d (e) is finite as e — > 0, then there exists 
a (possibly irregular) gauge transformation of the form 
(|2.8p such that the transformed one-parameter family of 
solutions satisfies assumptions 1, 2 and 3. 

Finally we note that, if one wanted to go to higher post- 
Newtonian orders, the assumptions used here to obtain a 
covariant approximation scheme would need to be mod- 
ified. First, to account for dissipative, radiative effects, 
one would need to introduce half-integer powers of e in 
the expansions. These would first arise at order 0(e 5 ^ 2 ) 
in g a b and order 0(e 7 / 2 ) in g ab . Alternatively, one could 
retain integer powers but make the replacement e — > e 2 
throughout. Second, it is well known that solutions of the 
post-l-Newtonian field equations are not good approxi- 
mations to exact solutions at distances > l/\/e. That 
is, although they work well in the near zone they break 
down in the local wave zone [22j. In order to find solu- 
tions which are good approximations everywhere, one has 
to match post-Newtonian solutions onto radiation zone 
post-Minkowskian solutions; see, for example, Blanchet 
[l|. However, the corresponding corrections to the near 
zone gravitational fields and to the dynamics of the bod- 
ies arises at post-2. 5-Newtonian order, and will therefore 
not be important for this paper. 

III. DERIVATION OF NEWTON-CARTAN AND 
POST-NEWTON-CARTAN THEORIES FROM 
GENERAL RELATIVITY 

In this section we derive the equations (|l.ip - f) 1 . 5[) 
of Newton-Cartan theory, (|1 .€>[) - (jl.lip of perturbative 
post-Newton-Cartan theory, and (|1.12al) - (I1.16P of com- 
bined post-Newton-Cartan theory, from the assumptions 
discussed in Sec. UU above. The derivation will be local 
and covariant 5 . 

We first note that it follows from assumption 1, to- 
gether with an adjustment of the normalization of the 



4 Dautcourt 19] shows that assumption 3 follows from assump- 
tions 1 and 2 when one makes additional assumptions about the 
global properties of the spacetime including asymptotic flatness. 
However, we will not follow this route here, since we want to 
obtain a purely local derivation of the Newton-Cartan and post- 
Newton-Cartan theories. 

5 A similar approach to deriving a covariant post-Newtonian the- 
ory was undertaken by L. Gunnarsen (unpublished) at the Uni- 
versity of Chicago in the 1980s. 
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one-form t a if necessary, that the covariant components 
of the metric can be expanded as 

9ab{£) = ~-t a t b + Pab + eqab + o(e), (3.1) 

where t a , p a b and q a b are e-independent tensor fields on 
spacetime. To see this, choose a basis of vector fields 
e% for a = 0,1,2,3 for which h ab = S rj e?et Then by 

1 3 

assumption we have A: 00 ^ 0, and in fact fc 00 must be 
negative in order for g a b(s) to have signature (— , +, +, +) 
for small e. Now expanding and inverting on this basis 
the expression (|2.2j) for the contravariant metric yields 
an expression of the form (|3.ip . 

Next, in any coordinate system we can compute the 
coefficients r^ 7 (e) of the connection by using the expan- 
sions (|3.ip and (|2.2[) of the covariant and contravariant 
metrics. This yields an expression of the form 

r^ 7 ( £ ) = Oie- 1 ) + O(e ) + 0(e) + o(e), (3.2) 

where the first three terms can be computed explicitly 
from the fields appearing in the metric expansions. Now, 
from assumption 3 it follows that the first term in Q3.2p 
vanishes, and the second term gives the coefficients of the 
Newtonian connection D a defined in Eq. (|2.5j) . Therefore 
we can write, for any one- form w a , 

V a (e)w b = D a w b - eA c ab w c + o(e), (3.3) 

where A c ab is an e-independent tensor field which is sym- 
metric in a and b. This quantity is the post-Newtonian 
perturbation to the connection. The e in brackets on the 
left hand side of Eq. (|3.3|) indicates the dependence of 
the derivative operator on e, not an application of the 
derivative operator to e. 

A. Newton-Cartan theory and perturbative 
post-Newton-Cartan theory 

1. Orthogonality conditions 

We start with the definition of the contravariant met- 
ric, 

g ab (e)g bc (e) = S a c , (3.4) 

and insert the expansions (|3.1[) and (|2.2p of the covariant 
and contravariant metrics. At order 0(e _1 ) this yields 
the Newton-Cartan orthogonality condition ([l.ip . and 
at order O(e ) it yields the corresponding post-Newton- 
Cartan condition (|1.6|) . 

2. Compatibility conditions 

Next, we insert the expansions (|3.ip and (|2.2p of the 
covariant and contravariant metrics and the expansion 
(I3.3P of the connection into the equations 

V a (s)g bc (e) = 0, V a {e)g bc {e) = 0. (3.5) 



At leading order, this yields the Newton-Cartan compat- 
ibility conditions (ll.2ap and (|1.2bp . and at subleading 
order one obtains the post-Newton-Cartan compatibility 
conditions (fj~T|) and ([Eg]) . 

3. Trautman conditions 

We next compute the Riemann tensor using the ex- 
pansion (|3.3I) of the connection. This yields 

Rabc^e(e)} = R abc d [D e ] +e[D b A d ac - D a A d bc ] 

+o(e). (3.6) 

Here the first term on the right hand side is the Rie- 
mann tensor of the connection D a , which we will denote 
henceforth simply as R ab d ■ From the symmetries of the 
Riemann tensor it follows that 

9 f[a m fM d] [Ve(e)}=0, (3.7) 

which is called the Trautman condition [l6[ 6 . Inserting 
the expansion (|2.2p of the contravariant metric and (|3.6p 
of the Riemann tensor, and expanding order by order 
in e, gives the Newton-Cartan Trautman condition (|1.3p 
at leading order, and the corresponding post-Newton- 
Cartan condition (ll.9p at subleading order. 

4- Stress-energy conservation 

Next, we insert the expansions Q3.3P and (|2.4p of the 
connection and stress-energy tensor into the conservation 
equation 

\7 a (e)T ab (e) = 0. (3.8) 

At leading order this yields the Newtonian stress-energy 
conservation equation (|1.5p . and at subleading order the 
post-Newtonian equation (jl.lip . 

5. Field equation 

Finally, we write the Einstein field equation (|2.ip in 
the form 

R a cb C [Ve(e)} - 47r[2g ac (e)g bd (e) - g ab (e)g cd (e)]T cd (e). 

(3.9) 

Inserting the expansions (|3.6p , (|3.1I) and (|2.4I) of the Rie- 
mann tensor, covariant metric and stress-energy tensor 
yields at leading order the Newton-Cartan field equation 
(|1.4p . and at subleading order the post-Newton-Cartan 
field equation (jl.101) . 



For many calculations it is advantageous to rewrite the Traut- 
man condition as g ea R eb ^ — 9 ed R ec t = 0. 
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B. Combined post-Newton-Cartan theory 

In this section we derive the equations (|1.12a[) - 
(|1.16p of combined post-Newton-Cartan theory from the 
Newton-Cartan and perturbative post-Newton-Cartan 
theories. 

Suppose we have a solution of the Newton-Cartan and 
perturbative post-Newton-Cartan theories, consisting of 
the fields t a , p ab , h ab , k ab , D a , A£ c , T ab and S ab . Such 
solutions posess a scaling symmetry corresponding to a 
change in units of time. Specifically, it is easy to check 
that for any real number A there is a mapping of solutions 
to solutions given by rescaling the fields by t a — > e x t a , 
k ab _^ e -2*k ab , T ab -> e- 4X T ab , S ab e- 6X S ab , Ag c 
e~ 2X A bc , with the other fields being left unchanged. If we 
now consider the expansions (|2.2[) . (|3.1[) . (|3.3|) and (I2.4[) 
of the contravariant metric, covariant metric, connection 
and stress-energy tensor, truncated to post-l-Newtonian 
order, we can apply this rescaling to effectively set e = 1 
in these expansions. Specifically, for any e > 0, we define 
the hatted fields by 

1 

ta = —pta-, (3.10a) 

h ab = h ab , (3.10b) 

Pab = Pab, (3.10c) 

k ab = ek ab , (3.10d) 
D a w b = D a w b - sA ab w c , (3.10e) 

fab = £ 2 T ab +£ 3 s ab^ (3 1Qf ) 

for any one- form w a . These will be the fundamental vari- 
ables of the combined theory. 

Next, we note that the orthogonality conditions (II. lj) 
and (|1.6|) are preserved under rescaling, which yields the 
orthogonality conditions (|1.12ap and (|1 . 12h>[) for the hat- 
ted variables. 

Next, using the definitions (|3.10j) of the hatted fields 
and the connection-compatibility conditions (II. 2[) . (II. 7[) 
and (|1.8[) we obtain 

D a (h bc + k bc ) = s 2 (A b ad k dc - A c ad fc bd ),(3.11a) 
Da(-i b ic+Pbc) = e(A d abPdc + A d acPbd ). (3.11b) 

The right hand sides of both of these equations are of 
post-2-Newtonian order. Therefore we can drop the 
right hand sides to obtain the compatibility conditions 
(I1.13[) ; this modification affects the theory only at post- 
2-Newtonian and higher orders, and not at Newtonian or 
post-l-Newtonian order. 

Finally, we can use exactly analogous arguments to de- 
rive the Trautman condition (|1.14l) , field equation (I1.15[) 
and stress-energy conservation equation (|1.16[) of the 
combined post-Newton-Cartan theory from the corre- 
sponding equations of the Newton-Cartan and pertur- 
bative post-Newton-Cartan theories. 



IV. DERIVATION OF STANDARD, 
COORDINATE-SPECIFIC POST-NEWTONIAN 
THEORY FROM PERTURBATIVE 
POST-NEWTON-CARTAN THEORY 

A. Change of viewpoint 

So far, we have shown that the Newton-Cartan and 
post-Newton-Cartan theories can be derived from general 
relativity together with the three assumptions discussed 
in Sec. M 

We now make a change of viewpoint, and consider 
these theories as independent theories in their own right, 
independent of general relativity. In other words, we for- 
get about the spacetime metric, and instead regard the 
fields t a , h ab , D a of the Newton-Cartan theory, and p ab , 
k ab and A c ab of the post-Newton-Cartan theory, as fun- 
damental. It is well known that the usual coordinate- 
dependent formulation of Newtonian gravity can be de- 
rived from the resulting Newton-Cartan theory, under 
the assumption of asymptotic flatness. In this section we 
will show that, similarly, the usual coordinate-dependent 
formulations of post-l-Newtonian theory can be derived 
from the post-Newton-Cartan theory, in suitably cho- 
sen coordinate systems, again under the assumption of 
asymptotic flatness. 

B. Derivation of Newtonian theory from 
Newton-Cartan theory 

We start by reviewing the well-known derivation at 
Newtonian order [l9j . We assume that the Newton- 
Cartan equations (jl.ll) - (jl.5[> are valid throughout all 
of spacetime, and that the Riemann tensor of the con- 
nection D a goes to zero at spatial infinity (asymptotic 
flatness). 

From the metric compatibility condition (ll.2bp . it fol- 
lows that there exists a function t on spacetime for which 
t a = D a t. We will call this function the time function. 
We now introduce a coordinate system x a = (x°, x 3 ) with 
x° = t. The orthogonality condition (jl.ip and Eq. (|1.2b[) 
then immediately lead to 

= = <5° (4-1) 

and 

= 0. (4.2) 

Next, the compatibility condition (jl.2bl) gives d^t v ~ 
r* v t\ = 0, where T^ v are the coefficients of the connec- 
tion D a , which yields 

C = (4.3) 

from Eq. (I4.ip . Similarly the compatibility condition 
(|I~2"aj) yields 

+ vi v h vP + K» hau = °- ( 4 - 4 ) 
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If we now define a spatial covariant metric h k j by 

h ik h kj = 8), (4.5) 
then Eq. (|4.4p results in 

T/ m = (hkLm + hkm.l — fym.fe) j (4.6) 

To; = 2^ ife (^w.o — e ki m B m ) , (4.7) 

and 

It = h* k <S> k . (4.8) 

Here the quantities $fc and i? m are still undetermined. 
Next, the field equations (II. 4[) imply 

i? l0 = 0, (4.9a) 
^ = 0, (4.9b) 

where i? a( g = R a p [D a ] is the Ricci tensor computed from 
the Newtonian connection D a . Computing the spatial 
components of this Ricci tensor explicitly, and simplify- 
ing using the condition (|4.3p . gives 

Ra = d k T% - diT k kj + rferk - r| fe r*. (4.io) 

Combining this with Eq. (|4.6p now gives Rij = 
( 3 > Rij[hki], where the right hand side is the three- 
dimensional Ricci tensor computed from the metric h k i 
at fixed t. Therefore, from Eq. (I4.9b|) . the Ricci tensor of 
the metric hij vanishes. Since this metric is three dimen- 
sional, it follows that the metric is flat. Hence at each 
fixed t we can choose the spatial coordinates so that 

hij=6ij. (4.11) 

Next, combining the field equation (|4.9a[) with the con- 
nection coefficients (|4.6[) - (14.8[) and simplifying using the 
coordinate condition (|4.11l) gives 

<, ,1,0,11, = 0. (4.12) 

Also the Trautman condition (|1.3[) implies that B is 
transverse, diBi = 0. Together with Eq. (|4.12j) this 
implies that the field B satisfies Laplace's equation, 
Bi,jj = 0. Next, using the assumption that R a p~ — > 
as r — > oo, we find that the only allowed nontriv- 
ial solutions to Laplace's equation are those with B = 
constant. These solutions can be eliminated by trans- 
forming to a uniformly rotating coordinate system. It 
follows that, in a suitably adjusted coordinate system, 
Bi = 0. Finally, the Trautman condition (|1.3[) implies 
that — = €ij k B k — 0, so that $i = <9i$ for some 
function $, which will be the Newtonian potential. It 
follows that 



while all the other connection coefficients vanish. The 
field equation (|1.4I) then reduces to the Poisson equation 
®,kk = 47rT 00 , and the stress-energy conservation equa- 
tion KT5J reduces to d a T aQ = 0, d a T m + $ ^T 00 = 0. 
These are the standard equations of Newtonian gravity. 

We note that the Newton-Cartan theory contains more 
local degrees of freedom that Newtonian theory. In par- 
ticular, if one assumes that Newton-Cartan theory holds 
only in a local region of spacetime, or if one assumes it 
holds everywhere but drops the assumption of asymptotic 
flatness, then one obtains a theory with an additional 
transverse vector field B that satisfies Laplace's equation. 
This is just the gravitomagnetic field which normally 
arises at post-l-Newtonian order. Thus, the Newton- 
Cartan theory admits source-free gravitomagnetic fields 
at Newtonian order. 

Another method that has been used in the literature 
to exclude these extra degrees of freedom is to assume 

that [n m, m 

h ec R abc d = 0, (4.14) 

where R abc d is the Riemann tensor of the connection 
D a . Augmenting the equations of Newton-Cartan the- 
ory with this assumption yields a covariant theory which 
is equivalent, locally, to Newtonian gravity. However, the 
assumption (|4.14p cannot be derived from general rela- 
tivity in a local manner; its validity requires the use of 
global information. For this reason we do not use the 
assumption (|4.14l) in this paper. 

C. Derivation of post-Newtonian theory from 
perturbative post-Newton-Cartan theory 

We now extend the above derivation to post-l- 
Newtonian order. We continue to use the adapted co- 
ordinate system derived above, and we assume that the 
post-Newton-Cartan equations are valid throughout all 
of spacetime. 

The spatial components of the orthogonality condition 
(|1.6[) imply that h lk p k j — <5*. Hence the covariant spa- 
tial metric h k j defined in Eq. (14.51) and the spatial com- 
ponents of p a b coincide, and in our adapted coordinate 
system we have 

Pij = h i:j = Sij. (4-15) 

The remaining components of the orthogonality condi- 
tion (HU) yield 

fc 00 = -1, (4.16) 

and 

k l ° = h lk p k0 = Pl0 . (4.17) 

In order to find expressions for pio and poo we have to 
consider the metric compatibility condition (jl.8[) , which 
yields 



(4.13) 



-poo, M + 2A MO p A0 . 



(4.18) 
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It follows that 

A° 00 = -^Poo,o + *iPio (4-19) 

and 

A% = -~Pooj- (4-20) 
The condition (jl.8l) also yields 

A° w = -Poi,„ + A^-pao + A^ oPiA , (4.21) 
which implies 

A%- = -po* (4.22) 

and 

A° 0j = $ j -p j0 , Q . (4-23) 

Since the connection is symmetric, Eq. (|4.22l) yields 
Poi,j — PQj.i = 0, and thus poi — g^, for some function 
g. Now under infinitesimal gauge transformations poi 
transforms as, from Eq. (IB4bl) . 

Poi -> Poi - ? ti , (4.24) 

where £° is an arbitrary function. Therefore by taking 
£° = g we can specialize the gauge to enforce 

PiH = 0. (4.25) 

Next, the symmetry of the connection applied to Eqs. 
P~2"0|) and flj^g]) yields together with Eq. (|4~2"5l that 
Poo,j = — 2$j, which implies poo = —2$ + x( x °)i where 
x{x ) is a function which depends only on x . Now note 
that the gauge condition (|4.25[) does not completely fix 
the gauge; from Eq. (I4.24[) gauge transformations with 
^° i = leave the condition (|4.25p invariant. Since poo 
transforms like poo ~ ^ Poo — 2£° under the infinitesimal 
gauge transformations of Eq. (IB4bl) , we can further spe- 
cialize the gauge to enforce 

Poo = -2$. (4.26) 

Simplifying Eqs. (|lT"9f . (R~2"0"|) and (R~2"2")) using the gauge 
conditions Eqs. (|4.25[) and (|4.26p now yields 

A° 00 = $ , A° i0 = $ i , A° y = 0. (4.27) 

Next we determine the remaining components of A c ab . 
From the compatibility condition (|1.7I) we get 

A^fcW + A j al h H = -*« Q - r^h* - (4.28) 
This leads to 

A% = -ifc« fc + W tfib (4.29) 

and 

A i oi =-i* y i0 + ^ (4-30) 



where Wyfe = —Wjik and = — Vji are undetermined. 
Since the connection is symmetric it follows from Eq. 
(!4~29l) that 

W m = W ikj + 1 (V% - fc* ,) . (4.31) 
Together with Wijk = —Wjik this implies that 

A< * = 5 (*V ( 4 - 32 ) 

Thus, from the metric compatibility conditions (| 1 . 8[) 
and (|1.7[) we have been able to determine all the compo- 
nents of A x ^ v , except for A l j and A l 00 . To determine 
those components we use the Trautman condition (|1.9p . 
which can be rewritten as 

k de R ecb a - k ae R ebc d = (D b A d ec -D e A d bc )h^ 
+ (D e A a bc — D c A a be ) h de . 

(4.33) 

Specializing to a = n, b = 0, c = k and d = m gives 

D A m n k - D n A m ok + D m A n ok - D c A\ m = 0, (4.34) 
which using Eqs. (|4.32p and (j4.30[) simplifies to 

V n k,m — V m k,n ~ Vnm,h = 0. (4.35) 

Since V nm is antisymmetric the solution of (|4.35[) is 

V mn = - (7 m ,„ - ln,m) , (4.36) 

where 7 m is some undetermined vector field. 

Next, considering the a — >■ n, b = c — > 0, d m 
components of Eq. (|4.33p . we find 

D A m n0 - D n A m oa + D m A n 00 - D A n 0m 

= \ (k ml Poo,m - k nl Poo , m i) ■ (4.37) 

If we define 

A m = A m 00 + l - Pw)tl k ml - 7m , , (4.38) 

then Eq. (14.37[) becomes 

A m , n - A n , m = 0, (4.39) 

so we can write 

An = ~l,m, (4.40) 

for some function 7. Combining Eqs. (|4.38[) and (|4.40[) 
now yields 

A ,; 00 = -^7,i + 7i,o - \pm,ik u , (4.41) 
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while Eqs. (j4~36|) and (jOU)) lead to 

^\ 3 = \{^-l^)-\k l \ - (4.42) 

We now define a metric g a b(s), for each e > 0, by the 
formula 

ffafc(ff) = —t a tb +Pab + Eq a b, (4.43) 
£ 



cf. the metric expansion (|3.1|) above. Here the tensor q a b 
is defined by 



9oo = 7, 
qoi = it, 
Qij = — k- 1 . 



(4.44a) 
(4.44b) 
(4.44c) 



We compute the inverse of this metric, using the values 
of the components of the fields t a , p a b and q a b given in 
Eqs. dUTJ, flUHJ) , 623 , (|4T2"d| and (|OI|) . The result is 
of the form [cf. Eq. (|2.2p above] 



sk 



ab 



£ 2 f b_ 



0(e 3 



(4.45) 



Here the components of the fields h ab and k ab (except for 
k ij ) are those given by Eqs. (O, (|4TTT|) . (jilB . gHH 
and (I4.25p . This result is guaranteed because we have im- 
posed the orthogonality conditions (jl.ip and (|1.6p . The 
fact that the spatial components of the coefficient of e 
in Eq. (|4.45[) are fc u follows from the choice (|4.44cj) of 
spatial components of q a b- The explicit form of the ten- 
sor j ab which appears in Eq. (|4.45[) will not be needed in 
what follows. 

Next, we compute the coefficients T"p (e) of the Levi- 
Civita connection associated with the metric (|4.43p , using 
the expansion (|4.45p . Suppressing indices, the result is 
schematically of the form 



f - -(htdt) + (hdp + ktdt) 

+e(hdq + kdp + jtdt) + 0(e 2 



(4.46) 



The leading order, 0(e~ r ) term vanishes identically by 
virtue of Eq. (14. ip . We can evaluate the next order, 0(e°) 
term using the the specific values of the components of t a , 
Pab, h ah and k ab in our adapted coordinate system, given 
by Eqs. gn}, BUI) . (|4~T51) . HPSt . (I4TT71) . (IT231) 

and (|4.26p . The resulting expressions are just the coef- 
ficients r^ 7 of the Newton-Cartan connection D a , with 
the only nonzero component being given by Eq. (|4.13p . 
Again, this result is not surprising because we have en- 
forced the compatibility conditions (|1.2a|) . (jl.2bp . (|1.7p 
and (|1.8p . Similarly, at the next order, we find that the 
0(e) components in Eq. (|4.46p coincide with the compo- 
nents of the post-Newton-Cartan field AJ3 , given by Eqs. 
(|4~27l) . ([4~32|) . (|OT|) and (|442)l . [Note that the result is 
independent of j , by virtue of Eq. (|4.ip ]. 

To summarize, we have been able to show that all of 
our Newton-Cartan and post-Newton-Cartan fields can 



be derived from the three fields t a , p a b and q a b that enter 
into the expansion (|4.43[i of the metric, using the stan- 
dard equations of general relativity. Moreover, that met- 
ric expansion is of the standard post-Newtonian form; 
using the specific values of the components of t a , p a b and 
qab given above and writing e = c -2 , Eq. (|4.43p takes the 
form 



ds 2 



2$ 
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O 



1 



dt 2 



+2 



dx l dt 



1 



dx l dx j . (4.47) 



This is the standard starting point for coordinate-specific 
post-Newtonian theory, involving a post-Newtonian cor- 
rection to the Newtonian potential $, and a gravito- 
magnetic potential % 7 . It follows that all of the re- 
lations of the Newton-Cartan and post-Newton-Cartan 
theories, when expressed in terms 7, 7* and gy, are 
either identically satisfied, or reduce to the Einstein equa- 
tions that one would compute directly from the metric 
(|4.43p . i.e. the coordinate-specific post-Newtonian equa- 
tions. Furthermore we have shown how to obtain the 
quantities $, 7, 7^ and qij, starting from solutions of 
the Newton-Cartan and post-Newton-Cartan equations. 
It follows that the equations of the perturbative post- 
Newton-Cartan theory are equivalent to those of the 
standard coordinate-specific post-Newtonian theory. 



V. DISCUSSION AND CONCLUSIONS 

We have derived a covariant version of the equations of 
the post-l-Newtonian approximation to general relativ- 
ity. These equations reduce to the standard coordinate 
formulation of post-Newtonian theory in asymptotically 
flat spacetimes in suitable coordinate systems. 

Although the covariant formulation is elegant, it does 
not provide a very compact or efficient representation of 
the theory. In a general coordinate system the combined 
post-Newton-Cartan theory involves 74 free functions to 
describe the geometry, as compared to four for stan- 
dard post-Newtonian theory. The covariant formulation 
is therefore mostly of formal interest. However, it may be 
useful to connect the different gauge-dependent formula- 
tions that are found in the literature. It may also be 
useful for deriving general properties of post-Newtonian 
theory. It might also provide insight into the meaning 
of the parameters of the parametrized post-Newtonian 



7 The spatial tensor qij does not contain any independent degrees 
of freedom; using the post-Newtonian field equations and making 



a gauge specialization gives qij = — 2<t>(5i 



This gauge spe- 



cialization is nearly always adopted in post-Newtonian theory. 
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framework [2|, if the analysis of this paper were gener- 
alized to the class of theories of gravity encompassed by 
that framework. 
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Appendix A: Limiting behavior of connection 
derived from other postulates 

In this Appendix we show that our assumption 3 on 
the limiting behavior of the connection will hold (up to 
gauge transformations) whenever the Riemann tensor is 
bounded as e — > 0. More precisely, whenever assumptions 
1 and 2 hold in a local region, and the Riemann tensor 
R abc d (e) is finite as e — > 0, then we show that there exists 
a (possibly irregular) gauge transformation of the form 
(|2.8p such that the transformed one-parameter family of 
solutions satisfies assumptions 1, 2 and 3 of Sec. Ill Al 

We start by fixing a coordinate system and computing 
the connection coefficients using the expansions (|2.2p and 
(|3.1|) of the contravariant and covariant metrics and the 
orthogonality condition (jl.ip . The result is of the form 
[cf. Eq. ((321) above] 

r^(s) = e-^- 1 ^ + r<°* + £ r«; 7 + 0(e 2 ), (Al) 
with 

r ( ~ 1} ; 7 - -/^ a (M[a^] + t p t [Xtil ) (A2) 

and 

r (0) ^ 7 = -&" A (M[A,/3] + *£*[A,7] + *A*(,8, 7 )) 

+ h aX (-pfa,x+Pp\,y+P7\,p)- (A3) 

Here r^ ^ are the coefficients of the Newtonian connec- 
tion D ai which were denoted simply f« in the body of 
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the paper. Also T ^ are the coefficients of the post- 
Newtonian connection perturbation, which were denoted 
A^ 7 in the body of the paper. We want to show that 

1 -^ 7 vanishes. 

We next compute the expansion of the Einstein tensor, 
which is of the form 

G^(e) = £ -2 G (-2)^ +e -l G (-l)a,3 +G (0)a / 3 



It follows from assumption 2 of Sec. Ill Al that the first 
four terms in this expansion all vanish. We find that 
G*(-2)a/3 vanishes identically, while G^~ 1 ^ Q ' 3 is given by 

G (-i) a p = 2H a ^H [i v h^ + ~h aP H' i u H v ll , (A5) 

where H a a = h^t^^y Setting this expression to zero 
yields E%W ^ = 0, from which it follows 19] that 



(A6) 



for some functions / and t. We now specialize the coordi- 
nates by choosing x° = t, so that h 0a = and t a = /<J°. 

We now extend this computation to the next order. 
The orthogonality relation (ll.6[) implies that 



fc 00 = -l// 2 , k 0i = h»p Qj /f, hv Pjk =5i. (A7) 



Using these relations and the expansions (|A2[) and (|A3[) 
gives G^ 00 = G^ Ql = and 

G (o)ij = G ij [hu] _ }_ D i D i f + -L/ ( ' •/)' /),,/•. (A8) 

Here the first term denotes the three-dimensional Ein- 
stein tensor computed from the metric hij = Pij (the 
inverse of h lJ ) 7 and Di is the covariant derivative asso- 
ciated with that metric. Also the O^ 1 ) piece of the 
Riemann tensor is given by 



R 



(-i) i 

OiO 



fDiD*f, 



(A9) 



with the other components being zero. Our assumption 
on the Riemann tensor forces this quantity to vanish, 
from which it follows from the vanishing of the expression 
(|A8|) for G^ 1 ^ that the Einstein tensor of the metric hij 
must be zero. We can therefore specialize the coordinates 
so that hij — — Sij . 

It now follows from Eq. (|A9I) that f^j = 0, so that 
/ = a(t) + Pi(t)x l . The leading order expression for the 
metric is therefore the Rindler metric, and we can ap- 
ply the standard gauge transformation 8 that takes the 
Rindler metric to the Minkowski metric. The result of 
this transformation is to effectively set / to unity, and so 



from Eq. ||£2) it follows that T ( ^ = 
formed one-parameter family of metrics. 



for the trans- 



In defining this gauge transformation we treat a and /3, as con- 
stants; their time dependence affects the final metric only at 
subleading order. The gauge transformation depends explicitly 
on e and is not smooth as e — y 0. 



Appendix B: Gauge freedom in the post-Newtonian 

fields 



In this Appendix we derive how the Newtonian and 
post-Newtonian fields transform under a regular gauge 
transformation ip E of the form (|2.9[) . Such a gauge trans- 
formation is parameterized by an e-independent diffeo- 
morphism ipo, and by a set of vector fields £m, £(2), • • -, 
one for each order in e. For simplicity, we will take tpo to 
be the identity mapping, since all quantities will trans- 
form trivially under this portion of the overall diffeomor- 
phism. Consider now any tensor field S(e) which depends 
on e, and has an expansion of the form 



(Bl) 



Here for brevity we have suppressed any tensor indices 
on S. We define the transformed expansion coefficients 
via the expansion 

tp E *S(e) = S{e) = S<°> + eg™ + £ 2 S (2) + 0(e 3 ). (B2) 

From Eq. (j!H)l) it now follows that [2l[ 

S<°> = S<®, (B3a) 
5(1) = S^+C^S^, (B3b) 
5(2) = + £? + C? 

+ i%%S<°>, (B3c) 

where C is the Lie derivative. 

We now apply this formalism to the expansions (12.21) . 
(|3.1[) , (|2.4p and (|2.5I) of the contravariant metric, co- 
variant metric, stress-energy tensor and connection. We 
use of the compatibility conditions (|1.2[) . denote gauge- 
transformed quantities with bars, and rewrite Lie deriva- 
tives in terms of D a derivatives. This yields that the 
Newtonian fields h ab , t a , D a and T ab are invariant, while 
the post-Newtonian fields transform as 

k ab = k ab -h ac D c £ b -h bc D4 a , (B4a) 

Pab = Pab - tJ c D b £ c - t b t c D a £, c , (B4b) 

A c ab = A c ab -2t; d R d{ab f + D {a D b) e, (B4c) 

S ab = S ab + i c D c T ab - 2T c{a D c £ b) . (B4d) 

Here we have written £m simply as £. One can check 
that the post-Newton-Cartan equations (|1.6[) - (|l.lip are 
invariant under these transformations, as they must be. 

To obtain the formula (|B4c[) . let uj b (e) be an arbitrary 
one form which depends smoothly on e, with the expan- 



sion UJ b 



,(0) 



eu b + 0{e 2 ). We define the tensor 



S a b = VaWf,, which has the expansion 



S ab (e) = S^+eS^+Oie 2 ) 



(B5) 



= D a ^ 



e D a u) : 



(i) 



A c a6 u;<°> +O&). 
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Applying the general transformation rule (|B3[) now yields 
SH^S^+C^, (B6) 

and 

(B7) 

Combining Eqs. (IB5|) - (|B7|) now yields 

+4 0) C a ^ (B8) 
which results in Eq. (|B4c[) . 

Appendix C: perfect fluids 

In this Appendix we describe as an example how per- 
fect fluids can be described in the covariant formalism. 
The stress-energy tensor is 

T ab = (p + p)u a u b +pg ab , (CI) 

where p is the density, p the pressure and u a the four- 
velocity. The appropriate form of the expansions of these 
fields is 

u a = Vi[< + £^ n + 0(e 2 )] , (C2a) 
p = ep a + e 2 p pn + 0(e 3 ), (C2b) 
p = £ 2 p n + £ 3 Ppn + 0(e 3 ). (C2c) 



Here the subscripts "n" and "pn" indicate the 
Newtonian-order and post-Newtonian order pieces. 
Comparing with the expansion (|2.4j) of the stress-energy 
tensor yields the formulae 

T ab = p n u a n u b n +pji ab , (G3a) 
S ab = (p pn + Pn )Kui + 2p n u^u%+ Pn k ab 

+P P nh ab . (C3b) 

Also the normalization of the four- velocity yields the con- 
ditions 

t a ul = 1, pab<ul = 1t a u a pn . (C4) 

One can check that inserting these expressions in the 
stress-energy conservation laws (11.51) and (ll.lip . using 
the specific forms of h ab , k ab and A£ c derived in Sec. 
IIVI and using the normalization constraints (|C4|) yields 
the usual equations of Newtonian and post-Newtonian 
hydrodynamics. 

Alternatively, one can combine the Newtonian and 
post-Newtonian pieces together, as in the combined post- 
Newton-Cartan theory derived in Sec. IIII B[ Defining 

P = Pn + Eppn, P = Pn + £Ppn and U a = < + £u£ n , then 

the combined stress energy of Sec. IIII Bl is 

t ab = (p + p)u a u b + p(h ab + k ab ), (C5) 

and the normalization constraint is u a u fc (t a ff,—p a fc) = — 1. 
Again one can check that these expressions lead to the 
usual post-Newtonian hydrodynamic equations. 



